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1 Abstract We investigate the minimal size of small superconducting grains by means of a 

t-H ' Ginzburg-Landau model confined to a sphere of radius R. This model is supposed to describe 

a material in the form of a ball, whose transition temperature when presented in bulk form, To, is 
' _h " known. We obtain an equation for the critical temperature as a function of 7? and of To, allowing 

us to arrive at the minimal radius of the sphere below which no superconducting transition exists. 
q '. PACS Number(s): ll.10.-z, 74.20.-z, 74.81.Bd 

U • 

5a ' I. INTRODUCTION 

S Recently, experimental results on small metallic grains [1] led to an important effort on theoretical investigations 
of superconducting pairing correlations in nanograins [2]. Many of the present developments in this area are based 
' on the exact solution of the discrete BCS model [3]. This model is described by a Hamiltonian for a discrete set of 
doubly-degenerate energy levels containing a pairing interaction for the scattering of pairs of electrons at levels next 
to each other. One of the fundamental questions addressed by these works has been stated long ago [4] and was 
O , restated in the following form [1]: What is the lower size limit for the existence of superconductivity in small grains? 

In the present letter we present an answer for the above question from another point of view, stemming from the field 
^-H [ theory framework of the Ginzburg-Landau model. Our approach consists in considering an Euclidean massive (\ip )d 
■ model describing a system constrained to be confined to a sphere of radius R. The rationale for our procedure is twofold, 
on one hand the Ginzburg-Landau model provides a well established and elegant theory of the phenomenology of the 
superconducting state. On the other hand, for Euclidean field theories, temperature, understood as imaginary time, 
and spatial coordinates share the same footing, which allows us to apply to the spatial coordinates a compactification 
mechanism using the Matsubara formalism [5-7]. We emphasize that we are considering an Euclidean field theory 
in D purely spatial dimensions, so we are not working in the framework of finite temperature field theory. Here, 
temperature is introduced in the mass term of the Hamiltonian (1.1) below by means of the usual Ginzburg-Landau 
prescription. 

We must stress here that we do not pretend to attain the degree of understanding coming from microscopic theories 
like BCS [8] or the reduced BCS [3], nevertheless we believe that our perspective can shed some light on this interesting 
1 question. 

We consider the Ginzburg-Landau Hamiltonian density, 



n= 1 -(V^ 2 + 1 -m 2 ^ + ^\ (1.1) 



in Euclidean 3-dimensional space, where u is the coupling constant, and m\ = a(T — To) is the squared bare mass (To 
being the bulk transition temperature of the superconductor and a > 0). From a physical point of view, introducing 
temperature by means of the mass term in this Hamiltonian, this should correspond to a spherical sample of material. 
We investigate the behaviour of the system as a function of the radius R of the confining sphere, in the approximation 
in which gauge fluctuations are neglected, and using spherical coordinates. Under these conditions one can write the 
generating functional of the correlation functions as, 



o 



Z = \ Vy*V(pexp \ dr driH{ip,Vip) , (1.2) 



with the field ip(r,8,(f>) satisfying the condition of confinement along the r-axis, ip(r = 0) = ip(r = R) = 0. These 
conditions of confinement of the r-dependence of the field to a sphere of radius R, permit us to proceed with respect 
to the r-coordinatc, in a manner analogous as it is done in the imaginary-time Matsubara formalism in field theory. 
The Feynman rules should be modified following the prescription, 



dk r 1 ^ , 2nn 

n— — oo 

The letter is organized as follows. In the next section we calculate the compactificd effective potential obtaining its 
R-depcndcnce. The suppression of superconductivity as a function of the grain radius is shown in section 3. Section 
4 present some concluding remarks. 

II. THE fl-DEPENDENT EFFECTIVE POTENTIAL 

We start from the expression for the one-loop contribution to the effective potential in absence of boundaries in 
spherical coordinates Ui(<po,R = oo), where in order to deal with dimensionless quantities in the required rcgular- 
ization procedures, we introduce parameters c 2 = m 2 / Air 2 a 2 , (i?/i) 2 = a -1 , g = (u/87r 2 ), (ifo/n) — <fto, ¥>o being the 
normalized vacuum expectation value of the field (the classical field) and [i a mass scale. In terms of these parameters, 



cm*, * _„)_^gfci£, # (2 .i) 



here m is the renormalized mass and the integration over the solid angle d f2 gives simply a factor 47r. Performing the 
replacements defined in (1.3), we obtain the boundary-dependent (i?-dependent) one-loop contribution to the effective 
potential in the form, 



8=1 
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This equation can be rewritten using the property 

+ °° 2 i J 

Y 7 J 1 = -T^Af(s-l;a), (2.3) 

(an 2 + c 2 ) s s-lda 1 y J K ' 
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where Af (s — 1; a) is one of the Epstein-Hurwitz zeta -functions [9] defined by, 

Ai ^ 01 ' ^ = ( fll n 2 + ... + arf n 2 + c 2 )^ (2 ' 4) 

valid for Re(v) > d/2 (in our case Re(s) > 1/2). Then Eq.(2.2) becomes, 

u^o, R ) = -4V« 3/2 E y^r/^ 8 x 

S— 1 ^ ^ 

xA A f( a -l;a). (2.5) 

The Epstein-Hurwitz zeta-function can be extended to the whole complex s -plane generalizing as is done in [5] the 
mode sum regularization procedure described in Ref. [9]. We write, 



1 — 1 TLi — 1 
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A$(u; ai ,...,a d ) = i; + 2^^(^ + c 2 r + 
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+2 d (am 2 1 + --- + a d n 2 d + c 2 y u . (2.6) 



Using the identity, 



we get, 
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dt f-V** 



(2.7) 



A° d (v;ai, ...,a d ) 



= j_r 
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dt t u ~ L e 
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where, 



ra» = l 

Tj(t,ai, ...,a,j) = Tj_i(i,oi, ...,aj_i)Ti(t,aj) , j = 2, 
Considering the property of functions Ti , 



ri(t,ai) = -\ + 



where 



(2.8) 

(2.9) 
(2.10) 

(2.11) 
(2.12) 
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we can notice that the surviving terms in Eq.(2.6) are proportional to (ai • • • a<j) t 1 / 2 ). Therefore we hnd, 
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A° d (i/;ai, ...,a d ) 
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i<j=l 



i=l 



(2.13) 

Inserting in Eq.(2.13) the explicit form of the function S(x) in Eq.(2.12) and using the following representation for 
Bessel functions of the third kind, K u , 

_ poo 

2(a/b)*K v (2Vab) = / dx x^e'^' 
Jo 

we obtain after some long but straightforward manipulations, 

Af(u;ai, ...,a d ) 
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(2.15) 



For d=l, taking v = s — 1, and identifying the compactified dimension with the radial coordinate, we obtain in our 
case from Eq.(2.5) the one-loop correction to the effective potential, 

(_!)«+! 2s 2 s - 1 /2 7r 2s-5/2 
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(2.16) 



or, remembering a 1 = (i?/z) 2 and </>o = Po/ft, 



U 1 (< Pa ,R) = J2 



.^+1/2^-3/2^2* 



n=l 
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where if_ s+ 3 are Bessel functions of the third kind. 



(2.17) 



III. CRITICAL BEHAVIOUR 



As we are analizing the case where the field has only one component we can neglect the Independence of the 
coupling constant, that is we consider u as the renormalized coupling constant. In this case, it is enough for us to use 
only one renormalization condition, 



-^-jtMvo, R)\ va =o 



m 



(3-1) 



Notice that we are using a modified minimal subtraction scheme [10], where the mass (and coupling constant, if it 
is the case) counter-terms are for even space dimension D, poles of Gamma-functions at the physical values of s 
(s = 1 for the mass, s — 2 for the coupling constant). For arbitrary D, we would have in the i?-independent term in 
Eq.(2.17), instead of the factor T(s — |), a factor T(s — -j), which would generate a pole at s = 1 for even dimensions. 
This polar term should be subtracted, giving the independent correction to the renormalized mass proportional to 
the regular part of the analytical extension of the Epstein-Hurwitz zeia-function in the neighbourhood of the pole at 
s = 1. For odd space dimensions, as for our case D = 3, there are no poles of Gamma-functions, but in order to have 
a coherent procedure in any dimension, we also subtract the corresponding term, performing a finite renormalization. 
Thus the independent renormalized mass, at one-loop approximation, is given by 



m 2 (i?) = ml + 



3\/2^ 3 / 2 
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(3.2) 



On the other hand, if we start in the ordered phase, the model exhibits spontaneous symmetry breaking, but for 
sufficiently small values of T" 1 and R the symmetry is restored. We can define the critical curve C(T C ,R) as the 
curve in the T x R plane for which the inverse squared correlation length , £~ 2 (T, R, ipo), vanishes in the independent 
gap equation, 



u f dil oj" n 



(3.3) 



where ipo is the normalized vacuum expectation value of the field (different from zero in the ordered phase). In the 
disordered phase, in particular in the neighborhood of the critical curve, ipo vanishes and the gap equation reduces to 
a i?-dependent Dyson-Schwingcr equation, 



m 2 (T,R)=m 2 (T) + 
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After steps analogous to those leading from Eq.(2.2) to Eq.(3.2), Eq.(3.4) can be written in the form 

m 2 (T,R) = ml 



(3.4) 
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Ki(nRm(T, R)) 



(3.5) 



If we limit ourselves to the neighborhood of criticality, m 2 (T, I) k 0, we may investigate the behavior of the system 
by using in Eq.(3.5) an asymptotic formula for small values of the argument of Bessel functions, 



K v {z)k\t{v)(^)~ V (*~0). 



Performing the derivative in Eq.(3.5) with the help of the formula, 

-^K v (z) = vz~ 1 K v (z) - K v+ i{z), 



(3.6) 



(3.7) 



and using the tabulated values T(l/2) = ^/tt, T(3/2) = y / 7r/2, we obtain after some manipulations in the neighbour- 
hood of criticality, 



(3i 



where C( z ) is the Riemann zeia-function defined for z > 1. 

Eq.(3.8) is of course meaningless as it stands. In order to obtain a critical curve and give a meaning to eq.(3.8), we 
perform an analytic continuation of the zeia-function ((z) to values of the argument z < 1, by means of its reflection 
property, 



(3.9) 



which defines a meromorphic function having only one simple pole at z = 1. The physical interpretation is achieved 
through a mass renormalization procedure that can be done as follows: remembering the formula 



lim 



1 



1 - z 



7, 



(3.10) 



where 7 = 0.57216 is the Euler constant, the renormalized mass m is then defined as 



m 2 (T,R) = lim 



\T,R) 



6nR(z - 1) 
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6^R' 



(3.11) 



Taking this renormalized mass equal to zero leads to the critical temperature, given as a function of the radius R by, 



(3.12) 



In Eq.(3.12), T corresponds to the transition temperature for the material in absence of boundaries (R — > oo), 
that is, to the bulk transition temperature. We see then that, in a spherical sample made of the same material, 
the critical temperature is diminished by a quantity proportional to the inverse of its radius. Also, we see that 
there is a minimal radius R^ below which superconductivity is suppressed. Identifying the usual 3-dimensional 
Ginzburg-Landau parameter (3 — u, the minimal radius is given by, 

R(0> = <& <« 3 » 



IV. CONCLUDING REMARKS 



In this letter, within the framework of the Ginzburg-Landau model, we were able to describe a spherical system, a 
grain, and we showed that, below a specific size, superconductivity is suppressed. Our method made use of effective 
potential and dimensional renormalization techniques supplemented by the Matsubara prescription (1.3) which give us 
a structure to work in a spherical geometry. We did not attempt to deal with the microscopic properties of ultrasmall 
metallic grains, however we believe that our approach can pave the way to further studies of the superconducting 
state in different geometries, notably through Monte Carlo and variational techniques. 
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